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ABSTRACT 

Multi-color light curves and radial velocities for TYC 1031 1262 1 have been ob- 
tained and analyzed. TYC 1031 1262 1 includesaCepheidwithaperiodof4.15270±0.00061 
days. The orbital period of the system is about 51. 2857±0. 0174 days. The pulsation 
period indicates a secular period increase with an amount of 2.46±0.54 min/yr. The 
observed B, V, and R magnitudes were cleaned for the intrinsic variations of the pri- 
mary star. The remaining light curves, consisting of eclipses and proximity effects, 
are obtained and analyzed for orbital parameters. The system consists of two evolved 
stars, F8II+G6II, with masses of Mi=1.640±0.151 M and M 2 =0.934±0.109 M 
and radii of Ri=26.9±0.9 R© and R2=15.0±0.7 R Q . The pulsating star is almost 
filling its corresponding Roche lobe which indicates the possibility of mass loss or 
transfer having taken place. We find an average distance of d=5070±250pc using the 
BVR and JHK magnitudes and also the V-band extinction. Kinematic properties and 
the distance to the galactic plane with an amount of 970 pc indicate that it belongs 
to the thick-disk population. Most of the observed and calculated parameters of the 
TYC 1031 1262 1 lead to a classification of an Anomalous Cepheid. 

Key words: stars: binaries: eclipsing - stars: fundamental parameters - stars: bina- 
ries: spectroscopic - stars: variables: Cepheids 
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1 INTRODUCTION 



Pulsating variables are located in a restricted region, so-called Cepheid instability strip (IS) in 
the Hetrzsprung-Russell (HR) diagram. This region includes 5 Scuti, RR Lyrae, Population I and 
II Cepheids which represent an important testing basis for theories on stellar structure and evolu- 
tion as well as pulsation mechanisms. These variables have very different masses, effective tem- 
peratures and chemical abundances. The stars in the IS are found not only in a wide range of 
masses but also found among very young Population I type stars to very old Population II type 
stars. These stars in the IS reveal many properties of the stellar interiors. 

At the beginning of the 20 th century Miss Leavitt made a great discovery that there is a relation 
between the period and a bsolute magn itude of the Cepheid variables in the Small Magellanic 
Cloud. Three decades later lBaadd (|1953|) called attention to the difference in the period-luminosity 
(hereafter P-L) relations of the Cepheids in the globular clusters and the classical Cepheids. Thus, 
the Cepheids are divided into two groups which are Population I (usually called Type I Cepheids) 
and Population II Cepheids (or Type II Cepheids), the formers are abou t 1 .5 mag brighter than the 
latter for the same pulsation periods. Taking into account this difference iBlaauw & Morganl (|1954|) 
made a revision on the P-L relation for the classical Cepheids. These variables are still used as 
distance indicators not only for our galaxy but also for nearby spiral and irregular galaxies. While 
the Type I Cepheids are young-disk stars associated with spiral arm and thin-disk populations 
evolved from main- sequence stars with masses from 3 to 15 solar masses, Type II Cep heids can 



be ol d-disk, thick-disk, or halo stars with masses somewhat smaller than the Sun (see 



Wallerstein 



(|2002l) . and references therein). Type II Cepheids are divided into three subclasses regarding their 



pulsation periods: BLHer stars with a peri od of < 7d, W Vir sta rs with 7< P < 20 d and RV 



Tau type variables with P> 20 d. Recently 



Soszvnski et al. 



(|2008|) suggested divisions at 4 and 



20 days. Masses and radii for these stars were estimated from their p ulsation properties. BL Her 



Harris 



(1981) showed that 



stars are common in globular clusters and in the disk of our galaxy, 
four of the 12 field BLHer stars have [Fe/H] values less than -1.5 which were identified as halo 
objects with a mean distance from the galactic plane of 1.8 kpc. In contrary, the 8 stars have 
metallicities near or above solar values. Their mean distance from the Galactic plane was about 
0.25 kpc, indicating characteristic of thick or thin-disk population. There are some metal-poor 
variables with periods between 0.5 and 3 or more days lying above the RR Lyrae-type stars which 
may have entirely different origin as compared to the normal Type II Cepheids in the globular 
clusters. These variables are called as Anomalous Cepheids (hereafter ACs). They are found in 
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the general field, globular clusters and nearby dwarf spheroidal galaxies. Anomalous Cepheids are 



brighter than the 



Soszynski et al. 



Type I 



2008) 



Cepheids at fixed colo r, and hence they follow a different P-L relation 



Fiorentino et al. 



d2006) (and references therein) suggested that the ACs 
are extension of the Type I Cepheids to lower metallicities and masses. Moreover, they predict a 
mass range for the ACs as 1.9<M/M <3, which may be a product of a mass about M< 4 M , 
experienced the helium flash in the past. 

Only six Type II Cepheids were known to be spectroscopic binaries. TXDel, IXCas, AUPeg 
and ST Pup are single-lined binaries for which spectroscopic observations and the resultant or- 
bits were published. Due to low orbital inclination no eclipses were detected in these binaries. 
Recently, two Ceph eid variables in eclipsing binaries , name ly T YC 1031 1262 1 a nd NSV 10993, 
were discovered by lAntipin. Sokolovski & Ignatieval (|2007l) and iKhruslovl (|2008l) . The pulsation 
periods are nearly identical with an amount of 4.2 d, but the or bital periods are about 5 1 and 40 
days. In the second part of the OGLE-III catalog of Variable stars ISoszynski et al.l(|2008|) presented 
197 Type II Cepheids of which seven are eclipsing variables and 83 ACs in the Large Magellanic 
Cloud. Since masses and radii for the stars were measured directly from the eclipse light curves 
and radial velocities, a detailed study of these Cepheids in the eclipsing binaries would be well 
worthwhile. 



The light variability of TYC 1031 1262 1 (ASAS J18261 1+12 
GSC 1031-01262, V=ll m 64, B-V=0 m .77) was discovered by 



2.6, 2MASS J18261 150+121234 9, 



Pojmanski. Pilecki & Syczygiel 



(12005b during the All-Sk y Automated Survey (ASAS-3). The va riability of this star was also 



detected independently by 



Antipin. Sokolovski & Ignatieval (|2007r) on Moscow archive plates in 



Crimea. They plotted the observations obtained by ASAS-3, photographic magnitudes obtained 
from the plates taken with the 40 cm astrog raph in the Crimea an d the observations obtained by 



the Northern Sky Variability Survey (NSVS, 



Wozniak et al. 



(|2004|) ) against the time and noticed a 



new variable as a Cepheid with some peculiarities. They started observations of TYC 1031 1262 1 
with the 50 cm telescope of the Crimean Laboratory. These CCD observations clearly revealed 
eclipsing nature of the star with a Type II Cepheid component in ou r Galaxy. The p ulsati on and 



Schmidt et al. 



(2009) made 



orbital periods are estimated to be 4.1523 and 51.38 d, respectively. 
VR photometry of Cepheid variable star candidates including TYC 1031 1262 1 in 2005 and 2006. 
They proposed a pulsation period of about 4.1508 d and classified it as a Type I Cepheid. 

In this study we present our multi-color photometric and spectroscopic observations of TYC 1031 1262 1. 
The main aim of this study is to derive the masses and radii of the components. Thus, mass and 
radius of a Type II Cepheid will be revealed, for the first time, directly from the spectroscopic and 
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photometric observations. We will also discuss the pulsation characteristics and its place in our 
Galaxy. 

2 OBSERVATIONS 

2.1 Photometric observations 

The photometric observations in the wide-band Johnson UBVR system were carried out with 
the 48 cm Cassegrain reflecting telescope and 35 cm MEADE LX200 GPS telescope at the Ege 
University Observatory. In the observations with the 48 cm telescope High-Speed Three-Channel 
Photometer and standard UBVR passbands were used. Thermoelectrically cooled ALTA U+42 
2048x2048 pixel CCD camera including BVR passbands was attached to Schmidt Cassegrain type 
Meade telescope. The BVR observations in 2008 were obtained on 63 nights between March 2 
and November 5. The observations in 2009 were obtained on 17 nights between August 3 and 
October 8. GSC 1031 193 and GSC 1031 1445 are taken as the comparison and check stars, re- 
spectively. Some basic parameters of the comparison stars taken from the SIMBAD database are 
listed in Table 1 . Although the program and comparison stars are very close in the sky, differential 
atmospheric extinction corrections were applied, especially for the observations obtained with the 
48 cm telescope. While all the program stars were observed simultaneously with the 35 cm tele- 
scope, each star was observed successively with the 48 cm telescope, i.e. in different times. The 
atmospheric extinction coefficients were obtained from observations of the comparison stars on 
each night. Moreover, the comparison stars were observed with the standard stars in their vicinity 
and reduced differential magnitudes, in the sense variable minus comparison, were transformed to 
the standard system. The standard stars are chosen from the lists of Landolt (1983, 1992). Helio- 
centric corrections were also applied to the times of the observations. 

In Fig. 1 we plot the U-, B-, V- and R-passband observations versus the pulsation period of 
the Cepheid variable. As it is seen the dominant light variations are originated from the pulsation. 
The below shifted observations correspond to the eclipses. The standard deviations of each data 
point are about 0.05, 0.03, 0.01, and 0.01 mag in U, B, V and R passbands, respectively. The 
observational data can be obtained from the authors. 

2.2 Spectroscopic observations 

Optical spectroscopic observations of TYC 1031 1262 1 were obtained with the Turkish Faint Ob- 
ject Spectrograph Camera (TFOSC) attached to the 1.5 m telescope on 6 nights (July-October, 
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Figure 1. The U-, B-, V-, and R-passband light curves for TYC 1031 1262 1. The ordinates and apsis are the standard magnitudes and the pulsation 
phases. 



Table 1. The coordinates, apparent visual magnitudes and the colors of the stars observed. 



Star 


a 


<5 


V (mag) 


B-V(mag) 


TYC 1031 1262 1 


\8 h 26 m 11.5 s 


12° 12' 34". 8 


11.64 


0.77 


GSC1031 193 


18 26 23.7 


12 15 47.2 


10.99 


0.66 


GSC1031 1445 


18 26 04.9 


12 07 11.6 


10.44 


0.63 



201 1) under good seeing conditions. Further details on the telescope and the spectrograph can be 



found at http://www.tug.tubitak.gov.tr. The wavelength coverage of each spectrum was 4000-9000 
A in 12 orders, with a resolving power of A/ A A 7 000 at 6563 A and an average signal-to-noise 
ratio (S/N) was ~120. We also obtained a high S/N spectrum o f the a Lyr (AO V) an d HD 50692 



2002) 



(GOV) for use as templates in derivation of the radial velocities (INidever et al. 

The electronic bias was removed from each image and we used the 'crreject' option for cosmic 
ray removal. Thus, the resulting spectra were largely cleaned from the cosmic rays. The echelle 
spectra were extract ed and wavelengths w ere calibrated by using Fe-Ar lamp source with help of 



the IRAF echelle (Tonrv & Davis 



1979) package. 

The stability of the instrument was checked by cross -correlating the spectra of the standard 
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star against each other using the FXCOR task in IRAK The standard deviation of the differ ences 



between the velocities measured using FXCOR and the velocities in iNidever et al. 
about 1.1 kms -1 . 



(120021) was 



3 PULSATION AND ORBITAL PERIODS 
3.1 Pulsation period 



he pu lsation period was estimated by Antipin. Sokolovski & Ignatieval (|2007l) and 



Schmidt et al. 



(120091) as 4.1523 and 4. 1508 d, respectively. Since the light variation of the pulsating component 



dominates in the light curve we fi rst attempted to refine the pulsation period using the program 



PERIOD04 (|Lenz & Bregeri 120051) . A Fourier power spectrum of all the available V-passband 
data gave a spectral peak at a frequency about / =0.2409 c/d which corresponds to 4.1581 d. 
This package computes amplitudes and phases of the dominant frequency as well as simultaneous 
multi-frequency sine wave fitting. Since the observations were obtained with different instruments 
in different years under dissimilar observing conditions, many systematic observational errors and 
computational errors affected the data. Therefore, we did not attempt to search for the second or 
the third frequencies for pulsation if they really do exist. As it will be explained in §4 we repre- 
sented all available V-data, subtracting the eclipses, with a truncated Fourier series which includes 
cosine and sine terms up to second order. We then calculated the light variation originated from the 
oscillations of the more luminous component. We separated all the data with an interval of about 
6 days and determined maximum times by shifting the calculated light curve along the time axis. It 
should be noted here that the shape of the light curve is assumed to be more or less constant during 
the time base of the observations. After the best fit is obtained, the times for the mid-maximum 
light are read off directly from the observations and presented in Table 2. Although the obser- 
vations obtained by ASAS and AAVSO have relatively large scatters we had to use all the data 
because of limited observations both in time elapsed and the c ontinuous observations due to its 



relativ ely longer pulsation period. Using the ephemeris given by 



(2007) 



Antipin. Sokolovski & Ignatieva 



Max(HJD) = 2 453 196.529 + 4 d .1523 x E (1) 

we obtained the residuals between the observed and calculated times of mid-maximum light 
as well as the number of the elapsed cycles. In Fig. 2 we plot the residuals O-C(I) versus the 
epoch numbers. The variation of the residuals resembles a parabolic change, in other words 
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Figure 2. The O-C residuals obtained by Eq.l for the pulsating star of TYC 1031 1262 1 are plotted against the pulsation cycle number. A parabolic 
fit to the data is shown b y solid line. The deviations from the parabolic fit are plotted in the bottom panel. Open circles refer to ASAS and AAVSO, 
and the dots to the data Antipi n. Sokolovski & I gnatieva ( 2007) and this study. 



TYC 1031 1262 1 appears to be undergoing gradual period increase. A least squares solution gives 
the following ephemeris, 

Max(HJD) = 2 453 196.201(0.061) + 4 d .15270(0.00061) x E + 9.73(2.13)10~ 6 x E 2 (2) 

The standard deviations in the last digits are given in the parentheses. In the bottom panel of 
Fig. 2 the deviations from the parabolic fit are also plotted. TYC 1031 1262 1 is undergoing period 
increase amounting to 2.46 (±0.54) min yr _1 . Despite the fact that the data base is very short, over 
five years, the characteristics of the variation in the pulsation period could be revealed. 



3.2 Orbital period 



The orbital period of TYC 1031 1262 1 was estimated by [Antipin. Sokolovski & Ignatiev a (2007) 
as 5 1 .38 d. We subtracted intrinsic variations of the more luminous star from all the available data, 
thus, the remaining light variations were assumed to be originated from the eclipses and proximity 
effects. We revealed light variations due to the eclipses and proximity analysis with an analysis of 
these data. The shape of the primary eclipse was revealed and compared by the observations which 
fall in the ascending and descending branches of the eclipse. Comparing the computed light curve 
with the observations, we obtained 28 times for mid-eclipse a nd presented in Table 3. The O 



C(I) residuals were computed using the first ephemeris given by I Antipin. Sokolovski & Ignatieva 
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Table 2. The times of mid-maximum light for TYC 1031 1262 1. The O-C(I) and O-C(II) residuals were 
computed with Eqs.l and 2, respectively. 



HJD-2 450 000 E O-C(I) O-C(II) Filter Ref. 



2710.533 


-117 


-0.177 


0.065 


V 


1 


2726.267 


-113 


-1.052 


-0.803 


V 


I 


2739.680 


-110 


-0.096 


0.158 


V 


1 


2785.286 


-99 


-0.166 


0.107 


V 


1 


2814.387 


-92 


-0.131 


0.152 


V 


1 


2859.573 


-81 


-0.619 


-0.323 


V 


1 


2885.053 


-75 


-0.053 


0.250 


V 


1 


2909.760 


-69 


-0.260 


0.049 


V 


1 


2909.773 


-69 


-0.247 


0.062 


V 


1 


3112.768 


-20 


-0.672 


-0.383 


V 


1 


3112.811 


-20 


-0.715 


-0.340 


V 


1 


3196.529 





0.000 


0.328 


V 


1 


3205.000 


2 


0.166 


0.494 


V 


2 


3213.529 


4 


0.391 


0.717 


V 


2 


3275.335 


19 


-0.088 


0.229 


V 


2 


3552.884 


86 


-0.743 


-0.521 


V 


1 


3553.291 


86 


-0.336 


-0.114 


V 


1 


3561.733 


88 


-0.198 


0.019 


V 


2 


3864.912 


161 


-0.137 


-0.126 


V 


2 


3868.925 


162 


-0.277 


-0.269 


V 


1 


3923.320 


175 


0.138 


0.099 


V 


1 


3939.820 


179 


0.029 


-0.026 


V 


2 


3943.926 


180 


-0.017 


-0.076 


V 


2 


3948.156 


181 


0.061 


-0.002 


V 


2 


4276.264 


260 


0.133 


-0.297 


V 


2 


4314.080 


269 


0.582 


0.099 


V 


3 


4342.949 


276 


0.385 


-0.138 


V 


3 


4355.414 


279 


0.393 


-0.148 


V 


3 


4388.986 


287 


0.747 


0.159 


V 


3 


4588.577 


335 


1.027 


0.130 


v 


4 


4588.627 


335 


1.077 


0.180 


R 


4 


4605.038 


339 


0.879 


-0.046 


R 


4 


4605.093 


339 


0.934 


0.009 


V 


4 


4638.476 


347 


1.099 


0.117 


V 


4 


4654.865 


351 


0.879 


-0.132 


R 


4 


4654.929 


351 


0.943 


-0.068 


V 


4 


4679.771 


357 


0.871 


-0.184 


R 


4 


4679.970 


357 


1.040 


0.015 


V 


4 


4684.179 


358 


1.127 


0.064 


V 


4 


4684.263 


358 


1.211 


0.148 


R 


4 



Ref: (1) ASAS, (2 ^Antipin. Sokolovski & Ignatieval J2007I) . (3) AAVSO (4) This study 



(I2007|) . In Fig. 3 we plot O-C(I) residuals versus the epoch numbers. A linear least squares fit gives 



the following ephemeris, 

Min(HJD) = 2 454 699.964(0.279) + 51 d .2857(0.0174) x E (3) 
The new orbital period is about 0.1 d shorter than the one estimated previously. 
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Table 3. The times of mid-minimum light for TYC 1031 1262 1. The O-C(I) and O-C(II) residuals were 
computed with the ephemeris given by Antipin, Sokolovski & Ignatieva 1 2007) and Eq. 3, respectively. 



rlj YJ-L hOU UUU 


e: 


U-t(l) 




Filter 


He J. 


3211.7872 


-29 


-8.520 


-0.892 


V 


1 


3212.8663 


-29 


-7.441 


0.187 


V 


1 


3213.3491 


-29 


-6.958 


0.670 


V 


1 


3211.6295 


-29 


-8.678 


-1.049 


V 


j 


3213.1512 


-29 


-7.156 


0.472 


V 


1 


3570.8060 


-22 


-6.562 


-0.873 


V 


1 


3572.1885 


-22 


-5.180 


0.510 


V 


1 


3571.5086 


-22 


-5.860 


-0.170 


V 


j 


3928.4599 


-15 


-5.969 


-2.218 


V 


1 


3929.3315 


-15 


-5.097 


-1.347 


V 


1 


3931.9234 


-15 


-2.505 


1.245 


V 


1 


3931.2237 


-15 


-3.205 


0.545 


V 


1 


3931.2279 


-15 


-3.201 


0.550 


V 


1 


3930.9966 


-15 


-3.432 


0.318 


V 


1 


3933.5099 


-15 


-0.919 


2.832 


V 


1 


4596.5385 


-2 


-1.003 


-0.854 


V 


2 


4598.4954 


-2 


0.954 


1.103 


V 


2 


4647.7061 


-1 


-0.844 


-0.972 


V 


2 


4647.7223 


-1 


-0.828 


-0.956 


V 


2 


4648.7162 


-1 


0.166 


0.038 


V 


2 


4649.0758 


-1 


0.525 


0.398 


V 


2 


4648.4924 


-1 


-0.058 


-0.186 


V 


2 


4650.1126 


-1 


1.562 


1.435 


V 


2 


4648.7959 


-1 


0.246 


0.118 


V 


2 


4699.6132 





0.054 


-0.350 


V 


2 


4699.6096 





0.051 


-0.354 


V 


2 


4750.8100 


1 


0.242 


-0.439 


V 


2 


4751.4869 


1 


0.919 


0.238 


V 


2 



Ref: m lAntipin, Sokolovski & Ignatievj |2007|) . (2) This study 



2. - 




-35 -30 -25 -20 -15 -10 -5 5 
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Figure 3. The O-C residuals obtained by Eq.3 for eclipsing binary TYC 1031 1262 1 and a linear lea s t squa res fit to the data. The deviations from 
the fit are plotted in the bottom panel. Open circles refer to the data Antipin, Sokolovski & Ignatieva (2007) and the dots to this study. 
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4 ANALYSIS 



4.1 Effective temperature of the primary star 



We have used our spectra to reveal the spectral type of the primary component of TYC 103 1 1262 1 . 
For this purpose we have degraded the spectral resolution from 7 000 to 3 000 by convolving 
them with a Gaussian kernel of the appropriate width, and we have also measured the equivalent 
widths (EW) of t he photospheric absorpt ion lines for the spectral classification. We have followed 
the procedures of iHernandez et al. I (I2004|) . choosing helium lines in the blue-wavelength region, 
where the contribution of the secondary component to the observed spectrum is almost negligible. 
From several spectra we measured -EW / Hei+FeiA4922 = 0.81 ± 0.07 A, £W H eI+FeIA4144 = 0.33 ± 
0.05 A, W Ca iA5589 = 0.35 ± 0.07 A, W Ca iA6i62 = 0.35 ± 0.07 Aand EWcma 



Hernandez et al. 



(20041), we 



1.05 ± 0.02 A. From the calibration relations £W-Spectral-type of 
have derived a spectral type of F8 II for the more luminous star with an uncertainty of about 1 
spectral subclass. In Fig. 4 we compare spectrum of the variable, obtained on JD2455 864, with 
the spectra of some standard stars. 

We have also observed the variable and comparison stars with the standard stars at the same 
nights. Peak to peak light and color variations of the pulsating star in V, B-V and V-R are about 
0.45, 0.14 and 0.09 mag, respectively. The average standard magnitudes and colors of the variable 
are obtained as <V> =11.64±0.04, <B-V>=0.77±0.06 and <V-R> =0.50±0.05 mag. Compar- 



ing the location of the variable in the (B-V)-(V-R) diagram given by Drilling & Landolt 



(2000) 



we estimate a spectral type of F8-9 II which is in a good agreement with that derived from spec- 
troscopy. The observed infrared colors of J-H=0.368±0.0 33 and H-K=0.10 6±0.033 are obtained 



using the JHK magnitudes given in the 2 MASS catalog (Cutri et al 



spond to a reddened F9±2 supergiant star (|Tokunaga 



2003). These colors corre- 



2000) which is consistent with that estimated 



both from the spectra and BVR photometry. 



Th e effective temp e rature deduced from the calib rati ons o f Drillin g & Landolt 



(119871) . 



Flower 



(119961) . 



Stra izvs & Kurielen e (1981) and 



Vitense ( 



(2000), 



de Jager & Nieuwenhu: 



1 98 1 ) are 5 750± 1 90 K, 5 740± 1 80 K, 



5 970±228 K, 6 000±150 K, and 5 830±180 K, respectively. The standard deviations are estimated 
from the spectral-type uncertainty. The weighted mean of the effective temperature was obtained 
for the primary star as 5 880±165 K. Th erefore, an interstell ar reddening of E(B-V)=0.21 mag is 



estimated from the tables given by lDrilling & Landolt 



(2000). 
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Figure 4. Comparison of the spectrum of TYC 1031 1262 1 with some standard stars with similar spectral type but different luminosity. 

4.2 Radial velocities 

To derive the radial velocities of the components, the 9 TFOSC spectra of the eclipsing binary were 
cross -correlated against the spectrum of GJ 182, a single-lined MOV star, on an order-by-order 
basis using the FXCOR package in IRAK The majority of the spectra showed two distinct cross- 
correlation peaks in the quadrature, one for each component of the binary. Thus, both peaks were 
fitted independently in the quadrature with a Gaussian profile to measure the velocity and errors 
of the individual components. If the two peaks appear blended, a double Gaussian was applied to 
the combined profile using the de-blend function in the task. For each of the 9 observations we 
then determined a weighted- average radial velocity for each star from all orders without significant 
contamination by telluric absorption features. Here we used as weights the inverse of the variance 
of the radial velocity measurements in each order, as reported by FXCOR. We adopted a two — 
Gaussian fit algorithm to resolve cross -correlation peaks near the first and second quadratures 
when spectral lines are visible separately. 

The heliocentric radial velocities for the primary (V p ) and the secondary (V s ) components are 
listed in Table 4 , along with the dates of observations and the corresponding orbital phases com- 
puted with the new ephemeris given in previous section. The radial velocities are plotted against 
the orbital phase in Fig. 5. The velocities in this table have been corrected for the heliocentric 
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Table 4. Heliocentric radial velocities of TYC 1031 1262 1. The columns give the heliocentric Julian date, the orbital phase (according to the 
ephemeris in Eq. 3), the radial velocities of the two components with the corresponding standard deviations. 



HJD 2400000+ Phase Star 1 Star 2 







v P 


(j 


V s 


a 


55747.4884 


0.4253 


4.1 


3.5 


29.7 


4.2 


55751.3090 


0.4998 


2.0 


3.3 






55796.2974 


0.3770 


-5.9 


5.4 


45.9 


5.6 


55800.3794 


0.4566 


5.4 


3.5 


22.7 


4.8 


55835.3002 


0.1375 


-7.7 


3.6 


48.8 


3.6 


55864.2074 


0.7011 


35.8 


3.4 


-35.9 


4.8 



reference system by adopting a radial velocity of 14 kms^for the template star GJ 182. The 
radial velocities listed in Table 4 are the weighted averages of the values obtained from the cross- 
correlation of orders #4, #5, #6 and #7 of the target spectra with the corresponding order of the 
standard star spectrum. The weight Wi = l/of has been given to each measurement. The standard 
errors of the weighted means have been calculated on t he basis of the e rrors (<Tj) in the velocity 
values for each order according to the usual formula (e.g. Topping (|l972)). The Oj values are com- 



Tonrv & Davis 



(119791) . The 



puted by FXCOR according to the fitted peak height, as described by 
observed radial velocities correspond to the pulsation phases of 0.718, 0.627, 0.360, 0.303, 0.608 
and 0.484. In these phases the variation in the radius of the primary star due to oscillation is small, 
therefore, radial velocity changes of the Cepheid caused by the pulsation are ignored. We did not 
attempt to decomposition the radial velocity measurements of the primary star into the pulsation 
radial velocity and the orbital radial velocity. 

First we analyzed the radial velocities for the initial orbital parameters. We used the or- 
bital period held fixed and computed the eccentricity of the orbit, systemic velocity and semi- 
amplitudes of the radial velocities. The results of the analysis are as follows: e=0.001 ±0.001, i.e. 
formally consistent with a circular orbit, 7= 1 1.07 ±0.85 km s -1 , i^i=27.4±1.7 and i^ 2 =48.1±1.7 
kms -1 . Using these values we estimate the projected orbital semi-major axis and mass ratio as: 
asim=76.50±2.44 R and q = f|=0.570±0.041. 



4.3 Intrinsic variations of the primary star 

The eclipsing binaries provide critical information about the orbital parameters such as orbital 
inclination, fractional radii, luminosities, ratio of effective temperatures etc. If the eclipsing binary 
is a double-lined binary, the masses and radii of the component stars can be determined in solar 
units. Using the inverse-square law one can accurately determine distance to the system which is 
independent of all distance methods. 
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Figure 5. Radial velocities folded on a period of 51.2857 days and the model. Dots (primary) and open circles (secondary) with error bars show the 
radial velocity measurements for the components of the system 



Table 5. Fourier coefficients of the oscillation light curves for TYC 1031 1262 1. 



Parameters 


B 


V 


R 


A 


1.0465±0.0022 


1.0316±0.0015 


1.0016±0.0012 


Ai 


0.2594±0.0028 


0.1948±0.0020 


0.1504±0.0016 


A 2 


0.0543±0.0028 


0.03S6±0.0020 


0.0271±0.0016 


Bi 


-0.0019±0.0032 


0.0123±0.0023 


0.0195±0.0018 


B 2 


-0.O530±O.0O31 


-0.0371±0.0022 


-0.0288±0.0018 



The observed light variations are composed of intrinsic light variations of the more luminous 
star and mutual eclipses. First we subtracted all the observations within the eclipses. The remaining 
observations are phased with respect to the pulsation period of 4. 15270 days. The light curves of 
the Cepheid show slow decline and rapid rise, i.e., asymmetric light curve with a sharp maximum. 
The median magnitudes are taken as 12.408, 11.639 and 11.144 mag for B, V and R passbands, 
respectively. The observed magnitudes are transformed to the flux using the median magnitudes. 
Then, we represented the intrinsic light variations of the primary star with a truncated Fourier 
series. Trial-and-error method showed that the observed light curves can well be represented by 
the second-order Fourier series. The coefficients are given in Table 5 and the fits are compared 
with the observations in Fig. 6. 

4.4 Analyses of the light curves 

The intrinsic light variations of the primary star were computed for each oscillation phase using 
the coefficients given in Table 5. After subtraction of the Cepheid light changes from the observa- 
tions we obtained light variations, consisting of only from the eclipses and proximity effects. In 
Fig. 7 the eclipsing light curves in the B, V and R passbands are plotted versus the orbital phases 
calculated by the ephemeris given in Eq. 3. The light curve of the system with curved maxima 
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Figure 6. The B, V, and R passbands light curves of the pulsating primary star and their Fourier representations (solid lines) with the coefficients 
given in Table 5. Note that the ordinates are normalized intensities. 

resembles to those (3 Lyrae-type binaries. The scatter in the binary light curves is mainly caused 
by the variations in the pulsation period. 

We have used the most rec ent version of the eclipsing binary light curve modeling algorithm 



of 



Wilson & Devinnev 



(|197 1|) (with updates, hereafter W-D), as implemented in the PHOEBE 
code of Prsa & Zwitter (2005). It uses the computed gravitational potential of each component to 
calculate the surface gravities and effective temperatures. The radiative characteristics of the stellar 
disks are determined using the theoretical Kurucz atmosphere models. The code needs some input 
parameters, which depend upon the physical properties of the component stars. The U-passband 
observations are limited and do not cover eclipses, therefore, we excluded these observations from 
the analysis of light curves for orbital solution. 

The BVR photometric observations for the system were analyzed individually. We fixed some 
parameters whose values were estimated from spectra, such as effective temperature of the hotter 
component and mass-ratio of the system, which are the key parameters in the W-D code. The 
effective temperature of the primary star has already been derived from various spectral type- 
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Table 6. Results of the BVR light curves analyses for TYC 1031 1262 1. 



v 



R 



Adopted 



"ff2 



«1 

S2 2 



^1 

X 
N 



74.85±0.13 
5880[Fix] 
5073±102 
3.416±0.037 
4.118±0.1S2 



71.26±0.39 
5880[Fix] 
4907±72 
3.594±0.026 
4.054±0.069 



76.31 ±0.49 
5880[Fix] 
4661±83 
3.709±0.038 
4.617±0.100 



74.0±0.4 
5880lFix] 
4840±80 
3.604±0.033 
4.292±0.104 



0.3615±0.0051 0.3380±0.0032 
0.1760±0.0101 0.2004±0.0049 
0.903±0.015 0.882±0.011 



2.532 
972 
0.0512 



2.186 
2557 
0.0293 



0.3258±0.0042 0.3378±0.0040 
0.1839±0.0060 0.1890±0.0064 
0.892±0.012 
0.653 
972 
0.0260 



effective temperature calibrations as 5 880 K and the mass-ratio from the semi-amplitudes of the 
radial velocity curves as 0.570. A preliminary estimate for the effective temperature of the cooler 
component is made using the depths of the eclipses in the BVR passbands. Therefore, initial linear 
and b olometric limb-dark ening coefficients for the primary and secondary components were taken 



from 



van Hamme ( 



1993|) . taking into account the effective tem peratures an d the wavelengths of 



the observations. The bolometric albedos were adopted from iLucy I (| 19671) as 0.5, typical for a 
fully convective stellar envelopes. The gravity-darkening exponents are assumed to be 0.32 for 
both components, because the stars are cool and are assumed to have convective envelopes. The 
rotational velocities of the components are assumed to be synchronous with the orbital one. We 
used Mode 2 of the W-D code which is for detached binaries with no constraints on the potentials. 
In this mode the fractional luminosity of the secondary component is computed from the other 
parameters via black body or stellar atmosphere assumption. We assumed synchronous rotation 
and zero eccentricity. 

The adjustable parameters in the differential correction calculation are the orbital inclination, 
the dimensionless surface potentials, the effective temperature of secondary, and the monochro- 
matic luminosity of the hotter star. Our final results are listed in Table 6 and the computed light 
curves (continuous line) are compared with the observations in Fig. 7. The secondary minimum ap- 
pears to be not very well reproduced. This is arisen mainly from the light variation of the pulsating 
star which is in front of the less massive secondary star at the secondary minimum. We assumed 
that light variation of the pulsating star is repeated with the same amplitude and the shape of its 
light curve is not change by time. The uncertainties assigned to the adjusted parameters are the 
internal errors provided directly by the Wilson-Devinney code. In the last three lines of Table 6 
the sums of the squares of residuals (x 2 ) , number of data points (N) and standard deviations (a) 
of the observed light curves are given, respectively. Taking weight inversely proportional to a, the 
adopted parameters given in the last column of Table 6 are obtained. 
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Figure 7. The B, V and R passbands light curves, originated only from the ecli pses and proximity effects, and the computed light curves. Circles 
and dots are for observations obtained by Antipin, Sokolovski & Ignatieva 1 2007) and in this study and solid lines are for the theoretical curves. 

5 RESULTS AND DISCUSSION 

A comparison of the results obtained in three light curves reveals some differences in the incli- 
nation angle and especially in the fractional radii of the stars. The B and R light curves are con- 
structed by only our own observations. The scatter in the B-passband observations is significantly 
larger compared to the other bands due to its faintness. Combining the spectroscopic results along 
with the photometric solutions, listed in the last column of Table 6, the absolute masses, radii, lu- 
minosities and surface gravities of the stars are obtained and presented in Table 7. The mass and 
radius of the pulsating component were determined for the first time directly from radial veloci- 
ties and multi-color light curves. They are unexpectedly large compared to the Type II Cepheids. 
The luminosity and absolute bolometric magnitude M^i of each star were computed from their 
effective temperatures and radii. The bolometric magnitude and effective temperature for the Sun 
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are taken as 4.74 mag and 5770 K ([Drilling & Landolt I (|2000|) ). Comparison of i ts location in th e 
Hertzsprung-Russell diagram (HR), constructed for the metal-poor cepheids by iGingoldl (| 19851) . 
their Fig. I, indicates that TYC 1031 1262 1 is in the instability strip of the Type II Cepheids, about 
20 times more luminous than the RR Lyrae stars . It seems to have the same luminosity and ef- 
fective temperature with the Type II Cepheids having pulsation periods of about 10 days, closer 
to the blue e d ge of the IS. When we compare with the solar composition models, for example 



Baraffe et al. 



1998), the primary star appears to be an evolved 5 M star and the secondary is 
consistent with an evolved 3 M star. 

After applying bolometric corrections to the bolometric magnitudes, we obtained absolute vi- 
sual magnitudes of the components. Taking into account the light contributions of the stars and 
total apparent visual magnitude we calculated their apparent visual magnitudes. The light contri- 
bution of the secondary star L 2 /(Li+L 2 )=0.02, 0.04 and 0.08 are obtained directly from the B-, 
V-, and R-bandpass light curve analysis, respectively. This result indicates that the light contri- 
bution of the less massive component is very small, indicating its effect on the color at out-of- 



eclipse is almost neg 



igible for the shorter wavelengths. For the bolometric corrections given by 



Drilling & Landolt I fo()Q(h ) we calculated the distance to TYC 1031 1262 1 as d=5074±207 pc. 



However the infrared JHK magnitudes give a distance of 5040±255 pc. Estimating distances to the 
stars are strongly dep ended on t he bolometric corrections. If w e adopt the bolometric corrections 



given by 



Code etaUl 1976) and 



Bessell. Castelli & Pled ( 



1998), the distances to TYC 1031 1262 1 



are obtained as 5082±540, and 51 10±650 pc, respectively. The average distance to the system was 
obtained as 5070±250 pc. The distance from the galactic plane is calculated as 970 pc using its 
galactic latitude of 11 degrees, about 14 times larger than the mean scale height of the Type I 
Cepheids. We have calculated Galactic space- velocity components of the system using the coordi- 
nates, distance, systemic radial velocity, and proper motions and listed in the last line of Table 7. 
The U, V, and W are velocities toward the galactic center, toward the direction of galactic rotation 
and toward the north galactic pole, respectively. With planer and vertical eccentricities of 0.50 and 
0.26 and the velocity components the system is mos t probably belonging t o the thick-disk popu- 
lation. This result confirms the earlier suggestion of lHarris & Wallersteinl (|1984|) that many field 
Type II Cepheids are metal rich and have kinematic properties similar to the old-disk or thick-disk 
population. 

The pulsation period of TYC 1031 1262 1 has shown a secu lar period increase with an amount 
to 2.46 (±0.54) min yr ^ rVinko. Szabados & Szatmaryl(| 19931) investigated secular period change 



in AU Peg, another Type II Cepheid which belongs to a binary system. The observed rate of the pe- 
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Table 7. Fundamental parameters of TYC 1031 1262 1. 



Parameter 



TYC 1031 1262 1 
Primary Secondary 



Spectral Type 
Mass (M Q ) 
Radius (R Q ) 
Teff (K) 
Luminosity (Lq) 
Gravity(cgs) 
a(R ) 
V 1 (kms" 1 ) 

1 

d (pc) 

/i a cos8, ^,5 (mas yr _1 ; 
U,V,W (kms' 1 ) 



F8(±1)II 
1.640±0.151 
26.9±0.9 
5 880±200 
764±144 
62±3 



G6(±1)II 
0.934±0.109 
15.0±0.7 
4 890±125 
109±26 
1 14±1 1 



79.58±2.54 
11.07±0.85 
74.0±0.4 
0.570±0.041 
5070±250 
2.7±0.8, 0.8±0.8 
-14.5±13.0, 47.7±14.9, -47.2±19.1 



riod in crease in TYC 1031 1262 1 is about 5 times larger than that of AU Peg. McCharty & Nemec 



1997|) observed a decrease in the pulsation period of V19 in the globular cluster NGC 5466 which 



is known as an AC with a mass of 1.66 M Q and a pulsation period of 0.82 d. They propose 
that this would be an indication of blueward evolution. In contrary, TYC 1031 1262 1 shows a 
period increase and its large magnitud e sho uld indica t e faste r evolutionary timescale redward. 



Sandage. Diethelm & Tammannl (| 19951) and 



Diethelm ( 



19961) proposed that all of the BL Her 



Cepheids appear to have increasing periods and are consistent with post-blue horizontal branch 
models of Population II stars. The observed period increase in TYC 1031 1262 1 appears to con- 
firm this suggestion. Unlike the BL Her Cepheids, the TYC 1031 1262 1 has a companion very 
close to it. Period change rates are determined for only a few ACs. Yet the cause of such a rel- 
atively large pulsation period change in a Type II Cepheid has not been obvious. Observations 
covering a long time span are crucially needed for better understanding of such a large increase 
rate in pulsation period in the TYC 1031 1262 1, as well as in other Type II Cepheids. 

As it is quoted in §1, Type II Cepheids are divided into three subclasses: BLHer, WVir, and 
RVTau. All three classes are characterized by the presence of Balmer emission, especially Ha, 
during some pa rts of their pulsation cycle. Evolutionary scheme of these classes is suggested by 
Gingoldl (|1985r) . In this scheme the BLHer stars are evolving from horizontal branch towards 
the lower asymptotic giant branch. However, the W Vir variables are on loops to the blue from 
the asymptotic giant branch. The RV Tau stars are moving to the blue in a post-asymptotic giant 
branch phase. He investigates evolutionary status of Type II Cepheids as post-horizontal-branch 
stars and calls attention to a different structure and evolution of the ACs. He estimates masses of 
Type II Cepheids around 0.6 with a radius of about 8 R© in the post-horizontal-branch phase. 
Bono. Caputo & Santolamazzal (|1997a|) suggest masses between 0.52 and 0.59 M Q for the Type II 
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Cepheids using their pulsation properties rather than evolutionary tracks. Now there is a consensus 
that Type II Cepheids are fundamental pulsators with masses below 0.8 M . 

With a mass of about 1.640 M Q , luminosity of 760 L and radius of 27 R Q the pulsating pri- 



mary component of TYC 1031 1262 1 cannot be a Type II Cepheid. '. 
mass would be greater than 4 M Q according to the models given by 



f it were a Type I Cepheid its 



Bono. Castellani & Marconi 



(12000 ) for Z=0.02 and even for Z=0.004. Howeve r the ACs are usually ac cepted to be metal-poor 



Soszvnski et al 



(120081) presented 83 ACs in the 



horizontal branch stars with masses above 1 .5 M Q . 
Large Magellanic Cloud and showed that ACs are located between Type I and Type II Cepheids in 
the P-L diagram. They are found in every dwarf spheroidal galaxy, and unlike Type I I Cepheids, 



are absent in globular clusters, except a few ACs in u Cen. iDemarque & Hirschfeldl (| 1975b pro- 
pose for their origin that they are y oung single stars due to recent star formation. In contrary, 
Renzini. Mengel & Sweigartl (|1977|) suggest that they are formed as a co nsequence of mass trans- 



Bono. Caputo & Santolamazza 



fer in b inary systems with the same age of the stellar systems they belong. 

1 1997b ) examined evolution and pulsation properties of ACs and concluded that models of masses 
between 1.5 and 2.2 M fit most of the stars. Concerning the origin of ACs, as single young 
stars or old binary systems, they have noted that their results did not support a single interpre- 
tation. The mass, luminosity, kinematic properties and light curves of TYC 1031 1262 1 fulfi l 



mos t of the properties of the A nomalous Cepheids. As pointed out by 



Soszvnski et al. 



and 



Fiorentino & Monelli 



(2008) 



(|2012r) the mode ide ntification of an AC is very complex and can not 



be based only on their light curves. However. ISoszynski et al. 



(|2008j) compared light curves of 



the ACs in the LMC. They proposed that the ACs pulsating in the first overtone have generally 
smoother light curves than for the fundamental-mode pulsators with rounded maxima and minima. 
The shape of the light curve of the pulsating star looks like to those of fundamental-mode p ulsators. 



Fiorentino & Monelli 



(120121) for the 



If we compare its location on the My - log P plane given by 
ACs in the LMC, whose pulsating periods are shorter than 2.4 d, the pulsating star is located on 
the extension of fundamental - pulsators, as though its classification as a fundamental-pulsator 
is supported. In Fig. 8 we c ompare position of TYC 1031 1262 1 on the c olor-magnitude diagram 



with other ACs taken from 



Harris. Olszewski & Wallerstein 



(1984) and 



Nemec. Nemec & Lutz 



(|1994|) . It is located among the brightest ACs close to the blue edge, with a luminosity of about 
800 times solar. The binary Cepheid ST Pup with a pulsation period of about 18.47 d, is located 
very close to the TYC 1031 1262 1. In contrary, the binary Cepheids TXDel with a pulsation pe- 
riod of 6.2 d and AUPeg with a period of 2.4 d, are located at the red-edge of the IS, with lower 
luminosity. 



© 0000 RAS, MNRAS 000, 000-000 



20 E. Sipahi et al. 



-5.1 
-4.1 
-3.1 

m 

J3 -2.1 
3 

-1.1 
-0.1 

0.9 

0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
(B-V)„ 

Figure 8. The position of TYC 1031 1262 1 on the color-magnitude diagram. The circles refer to the Type II Cepheids in globular clusters and the 
dots for TYC 1031 1262 1, ST Pup, AU Peg and TX Del. The lines show the borders of instability strip for Type II Cepheids. 

6 CONCLUSIONS 

We have obtained multi-color light curves and radial velocities of the eclipsing binary system 
TYC 1031 1262 1 of which the luminous component is a Cepheid. The astrophysical parameters 
of the component stars were obtained directly by analysing the light curves and radial velocities. 
In earlier studies the brighter star was classified as a metal-poor Type II Cepheids. In contrary, 
some investigators include it into the classical cepheids. The mass of the pulsating primary star 
with an amount of 1.64 M Q indicates that its structure and evolution are different from those of 
the RR Lyrae and W Vir stars. The galactic velocity components U, V, W have been determined as 
-14.5, 47.7 and -47.2 km s _1 and the distance from the galactic plane as about 970 pc. All these 
properties and the asymmetric pulsating light curve lead us to classify the Cepheid as an Anoma- 
lous Cepheid. The pulsation period appears to be increasing with an amount of 2.46 (±0.54) min 
yr~\ the largest period change observed in the Type II Cepheids so far. This observed increase rate 
indicates that the star evolves redward in the IS with shorter timescales. 

The pulsating primary star fills almost 85 per cent of its corresponding Roche lobe while the 
less massive star 61 per cent. The pulsating component in TYC 1031 1262 1 is very close to its 
Roche lobe, therefore, mass loss or transfer to its companion has taken place. Both components 
are supergiant and the separation is only two times the sum of their radii. The companion should 
unexpectedly affect both pulsation and evolution of the Cepheid. However, it is not clear how such 
a close companion can cause any period variation on pulsation by tidal interaction alone. Most of 
the ACs are faint stars and have relatively longer orbital periods, therefore, measurements of the 
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radial velocities ACs in binary systems are very difficult. In addition the velocity variations due 
to the orbital motion may be affected by the pulsation depending on the inclination of the orbits. 
In the case of TYC 1031 1262 1 the orbital inclination is sufficiently large to separate the orbital 
velocities of both components. If mass loss or transfer from the pulsating star has taken place, the 
orbital period of the system would be changed, which can be revealed by the observations to be 
made in coming years. This will yield us some clues about the further evolution of the ACs in 
binary systems. 
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